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SOLUTIONS TO THREE-DIMENSIONAL NAVIER-STOKES 
EQUATIONS FOR INCOMPRESSIBLE FLUIDS 

JORMA JORMAKKA 



' Abstract. This article gives explicit solutions to the space-periodic Navier- 

Stokes problem with non-periodic pressure. These type of solutions are not 
unique and by using such solutions one can construct a periodic, smooth, 
divergence-free initial vector field allowing a space-periodic and time-bounded 
external force such that there exists a smooth solution to the 3-dimensional 
Navier-Stokes equations for incompressible fluid with those initial conditions, 
. but the solution cannot be continued to the whole space. 

^' 



1. Introduction 

Let X = {xi,X2,X3) e denote the position, i > be the time, p{x,t) G K be 



e2 



the pressure and u{x,t) = {ui{x,t))i<i<s G be the velocity vector. Let fi{x,t) 
• be the external force. The Navier-Stokes equations for incompressible fluids filling 

ID '. all of for t > are [1] 

^ ■ 3 

^+^u,gi = «Azi,-^+/,(x,i), x6M^ i>0, l<z<3 (1.1) 

^ : 3 

div?/ = y^=0, xeR^t>0 (1.2) 
^— ^ ' oxi 

i=i * 

with initial conditions 

u{x,0) = u°{x), xeR3. (1.3) 
I Here A — J2'i=i the Laplacian in the space variables, w is a positive coefficient 

and vP{x) is C°°(M'^) vector field on required to be divergence-free; i.e., satisfying 
divu'^ = 0. The time derivative at t = in (1.1) is taken to mean the limit 
when t — J> 0+. 

This article shows that there exists C°°(M^), periodic, divergence- free initial 
vector fields vP defined at K.^ such that there exists a family of smooth (here, in 
the class C°°(R'^ x [0,cx)))) functions u{x,t) and p{x,t) satisfying (1.1), (1.2) and 
(1.3). We also show that there exist a periodic and bounded external force fi{x,t) 
such that the solution cannot be continued to the whole M"^ x [0, oo). 
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2. Theorems and Lemmas 

The simple explicit case of u{x,t) in Lemma 2.1 satisfies the conditions given 
in (1.3) and allows a free function g{t) that only satisfies g{0) = g'{0) = 0. The 
solution is then not unique. If the time derivatives of u{x, t) are specified at t = 
then the solution in the lemma is unique. 

Lemma 2.1. Let 

Ui = 27rsin(27ra;2) + 27rcos(27ra;3), 

U2 = 27rsin(27ra;3) + 27rcos(27ra;i), 

W3 = 27rsin(27ra;i) + 27rcos(27ra;2) 

be the initial vector field, and let fi{x,t) be chosen identically zero for 1 < i < 3. 
Let g -.M. M. be a sm,ooth function with g{0) = g'{0) = and /3 = (27r)^w. The 
following family of functions u and p satisfy (1.1)-(1.3); 

ui = e-^'2TT (sin(27r(a;2 + g{t))) + cos{2^[xz + g{t)))) - g'{t), 

U2 = e-f^'2TT (sin(27r(x3 + g{t))) + cos(27r(a;i + g{t)))) - g\t), 

Us = e-'^*27: (sin(27r(x'i + git))) + cos(27r(a;2 + g{t)))) - g'{t), 

p = -e-2^*(27r)2 sin(27r(a;i + g{t))) cos(27r(a;2 + g{t))) (2.1) 

- e-2'3*(27r)2 sin(27r(a;2 + g{t))) cos(27r(a;3 + g{t))) 

3 

- e"^'^*(27r)2 sin(27r(a;3 + g{t))) cos(27r(a;i + g{t))) + g"{t)'^Xj. 

Proof. The initial vector field is smooth, periodic, bounded and divergence-free. 

Let (i, k, to) be any of the permutations (1, 2, 3), (2, 3, 1) or (3, 1, 2). We can write 
all definitions in (2.1) shorter as (here g'{t) = dg/dt): 

= e-'^*27r (sin(27r(xfe + g{t))) + cos(27r(x„ + g{t)))) - g'{t), 

p = -e"2'^*(27r)2 sin(27r(xj + g{t))) cos(27r(xfc + g{t))) 

- e-2'3*(27r)2 sin(27r(a;fe + git))) cos(27r(.T„ + g{t))) 

3 

- e-2/3*(27r)2 sin(27r(a;„ + g{t))) cos(27r(a;i + g{t))) + g"{t) ^ Xj. 

It is sufficient to proof the claim for these permutations. The permutations (1, 3, 2), 

(2, 1. 3) and (3, 2. 1) only interchange the indices k and to. The functions (2.1) are 
smooth and u{x,t) in (2.1) satisfies (1.2) and (1.3) for the initial vector field in 
Lemma 2.1. We will verify (1.1) by directly computing: 

du ' 

= -j3e-^'2Ti {sm{2Ti{xk + g{t))) + cos(27r(a;„ + g{t)))) - g'\t) 
+ g'(t)e"'3*(27r)2 (cos(2^(a;fe + git))) - sin(27r(a;™ + g{t)))) , 
-vAui = ve~^\2'K)^ (sin(27r(a;fc + g{t))) + cos(27r(a;™ + git)))) , 

^ = -e-2/3t(27r)3 cos(27r(a;i + git))) cos(27r(a;fe + g{t))) 

+ e-2'3*(27r)3 sin(27r(a;„ + git))) sin(27r(a;i + git))) + g"it). 
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The functions Uk and Um are 

Uk = e-^'2TT (sin(27r(a;™ + g{t))) + cos(2^(.t, + g{t)))) - g'{t), 
Um = e"^*27r (sin(27r(xi + g{t))) + cos{2Tr{xk + g{t)))) - g'{t). 

The remaining term to be computed in (1.1) is 

Edui dui dui dui 



. ,. , , dxj dxi dxk dxm 

= e-'^^\2-Kf cos{2-K{xi + g{t))) cos{2-K{xk + g{t))) 

- e-^^*{2nf sm{2n{xm + g{t))) sin(27r(a;i + g{t))) 

- g'{t)e-^\2nf (cos(27r(a;fc + g{t))) - sm{2n{xm + git)))) • 
Inserting the parts to (1.1) shows that 

This completes the proof . □ 

thml Theorem 2.2. There exists a periodic, C°°(R'^), and divergence-free vector field 
vP{x) = {u^{x))i<i<3 on M."^ such that the following two claims hold: 

CI: The solution to (1.1)-(1.3) is not necessarily unique. In fact, there are 
infinitely many C°° {M? x [0,oo)) functions u{x,t) = {ui{x,t))i<i<s and 
p{x,t) satisfying (1.1), (1.2) and (1.3). 
C2: Periodic initial values do not guarantee that the solution is bounded. In- 
deed, there exist unbounded u{x , t) and p{x,t) satisfying (1.1), (1-2) and the 
initial values (1.3). There also exist bounded solutions that are periodic as 
functions of x. 

Proof. Let fi{x,t) be chosen identically zero for 1 < i < 3, and let us select g{t) = 
^ct^ in Lemma 2.1. The value c € R can be freely chosen. This shows CI. If 
c = then the solution is bounded and it is periodic as a function of a;. If c ^ 
then Ui{x,t) for every i and p{x,t) are all unbounded. In Ui{x,t) the ct = g'{t) 
term and in p the term c{xi + X2 + Xs) = g"{t) Y^^=i ^^"^ '^ot bounded. This 
shows C2. The failure of uniqueness is caused by the fact that (1.1)-(1.3) do not 
determine the limits of the higher time derivatives of u{x,t) when t 0+. These 
derivatives can be computed by differentiating (2.1) but the function g{t) is needed 
and it determines the higher time derivatives. As g{t) can be freely chosen, the 
solutions are not unique. □ 

thiii2 Theorem 2.3. There exists a smooth, divergence-free vector field u° {x) onM.^ and 
a smooth f{x, t) on x [0, oo) and a number Ca,m,K > satisfying 

M°(a; + ej) = u°(a;), f{x + ej,t) = f{x,t), l<j<3 (2.2) 

(here ej is the unit vector), and 

\d:drfix,t)\<CcmK{i + \t\r'' (2.3) 

for any a, m and K, such that there exists a > and a solution {p,u) of (1.1), 
(1.2), (1.3) satisfying 

u{x, t) = u{x + ej,t) (2.4) [eT] 
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on X [0, a) for 1 < j < 3, and 

P,ugC°°{R^ x[0,a)) (2.5) \78_ 

that cannot be smoothly continued to M.^ x [0, oo). 

Proof. Let us make a small modification to the solution in Lemma 2.1. In Lemma 
2.1, 5 : M — > M is a smooth function with g{0) = g'{0) = 0, but we select 

g{t) = -ct^ 

' 2 a-t 

where c ^ and a > 0. 

The initial vector field u'^(x) in Lemma 2.1 is smooth, periodic and divergence- 
free. The period is scaled to one in (2.1). The f{x, t) is zero and therefore is periodic 
in space variables with the period as one. Thus, (2.2) holds. The constant Ca,m,K 
is selected after the numbers a,m,K are selected. The force f{x,t) is identically 
zero, thus (2.3) holds. The solution (2.1) in Lemma 2.1 is periodic in space variables 
with the period as one. Thus (2.4) holds. The solution u{x,t) in (2.1) is smooth if 
t < a, thus (2.5) holds. The function g'{t) has a singularity at a finite value t = a 
and g'{t) becomes infinite a,t t = a. From (2.1) it follows that if t approaches a 
from either side, there is no limit to the the oscillating sine and cosine term in ui, 
and the g'{t) additive term approaches infinity. Thus, the solution u{x,t) cannot 
be continued to the whole M.'^ x [0, oo). This completes the proof. □ 

thmS Theorem 2.4. There exists a smooth, divergence-free vector field vP {x) onM.^ and 
a smooth f{x, t) on x [0, oo) defined as a feedback control function using the 
values of u{x,t) and a number Ca.m,K > satisfying 

u°{x + Cj) = u°{x), f{x + ej,t) = f{x,t), l<i<3. (2.6) 

(here Cj is the unit vector), and 

\dyrf{x,t)\<c„^K(i + \t\)-'' (2.7) ^ 

for any a, m and K, such that there exist no solutions {p,u) of (1.1), (1-2), (1.3) 
on X [0, oo) satisfying 

u{x,t) = u{x + ej,t) (2.8) [eT* 

on M.^ X [0, oo) for 1 < j <3, and 

e C°°(R-^ X [0,oo)). (2.9) 

Proof. Let the solution in Theorem 2.3 with the particular g be denoted by U and 
a be larger than 1. A feedback control force f{x,t) is defined by using the values 
of the function u{x, t') for t' < t. In practise there is a control delay and t' < t but 
we allow zero control delay and select f{x,t) as 

fi{x,t) = ^Ui{x,t) - ^Ui{x,t). 

Inserting this force to (1.1) yields a diff'erential equation in space variables 

dUi dui ^ dp T,,^ „ , 

at oxj oxi 

This force is defined in the open interval t G [0, a) and can be smoothly continued 
as zero to [a, oo). There is a solution u{x,t) = U{x,t) to this equation. We notice 
that if u{x,t) = U{x,t) then the force fi{x,t) takes zero value at every point. This 
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is correct: if wc apply external control force without any control delay, it is possible 
to keep the; solution u exactly at the selected solution U. This is not the same as 
to say that there is no force. If the solution u{x,t) would be different from U{x,t), 
then the force would not be zero. Since ui{x,t) = Ui{x,t) becomes infinite when t 
approaches a, the solution cannot be continued to the whole space x [0, oo). As 
in Theorem 2.3 the conditions (2.6)-(2.8) hold, but (2.9) does not hold. 

We must still discuss if the feedback control force can control the equation (1.1) 
or if there can be several solutions. When the higher order time derivatives of 
u{x, t) are fixed at t = the solution to (1.1) is unique because of the local-in-time 
existence and uniqueness theorem. This means that if a solution u{x, t) starts as 
U{x,t) in some small interval t < e for some small e < 1, then it will continue as 
U (x, t) for all times t < a ii the external force is zero. The question is whether the 
feedback control force f{x,t) can steer the solution u{x,t) to a possible solution 
U{x,t). The external force can freely change the higher order time derivatives of 
u{x,t) in the equation (1.1). Thus, the external control force can set the higher 
order time derivatives of u{x,t) to those of U{x,t), therefore the answer is that 
the external control force can control the equation and the solution can be set to 
U{x,t). 

The difference between Theorems 2.3 and 2.4 is that if the force f{x,t) is zero 
in Theorem 2.3, there is a family of solutions corresponding to different selections 
of g{t), but if the force f{x,t) is zero in Theorem 2.4, then necessarily the solution 
u{x,t) equals U{x,t) because otherwise the force is not zero. 

Let us mention that we may select a force that does not take the value zero at 
all points e.g. by adding a control delay that has a zero value at t = and when 
t > ti for some fixed ti satisfying < ti < a and smoothing the force to C°°. At 
some points t < ti the control delay is selected as nonzero and consequently the 
force is not zero at all points. This completes the proof. □ 



Let us continue by partially solving (1.1)-(1.3). Firstly, it is good to eliminate 
p by integrability conditions as in Lemma 2.5. We introduce new unknowns hi^k- 
The relation of hi^k and Ui is given by Lemma 2.6. 

Iem2 Lemma 2.5. Letu{x,t) andp{x,t) be C°° (R^ x [0, cxd)) functions satisfying (1.1) 
and (1.2) with fi{x,t) being identically zero for 1 < i < 3, and let {i,k,m) be a 
permutation of (1,2,3). The functions 

_dui duk 



satisfy 



'hi,k ^ ^ 0/t»,fc ^ _ dum ^ ^ dUm duj ^ dUm dUk (2 11) 

dt 4^ ^ dXj dXm dXk dXm dXi dXm 



for all X and i > 0. 

Proof As peC^ (M^ X [0, oo)), 

d dp d dp 
dxi dxk dxk dxi 
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Thus, from (1.1) when fi{x,t) is identically zero for 1 < i < 3, we obtain 

dp dui r d'^Ui duj dui \ ^ dui 

^\ dx^dxh dxhdxn/ < 



ddtdxk V dxjdxk dxkdxj/ dxk 

j=i ■' ■' 

_ dp_du^ / d^uk ^ duj_duk\ _ ^^duk 

ddt dxi ^\ ^ dxndxi dxidxj) dxi 



This yields 



3 3 3 

J, , A I, sr^dujduk duj dui 



ddt ' ^ dxj ' ^ dxi dxj ^ dxk dxj 



The right-hand side of (2.12) can be written in the form 

dui duk duk duk du,n duk du-i dui duk dui dum dui 



dxi dxi dxi dxk dx^ dxm dxi, dx^ dxk dxk dxk dxm 
dui / duk _ dui \ ^ duk / duk _ diH \ ^ dum dtn _^ dum duk 
dxi V dxi dxk ' dxk ^ dxi dxk ' dxk dxm dxi dxm 



In (2.13) we have replaced the sum X]j=i by X]je{2 km} which is possible since 
{i, k, m) is a permutation of (1, 2, 3). Due to (1.2) the right-hand side of (2.13) can 
be further transformed into 

dura ( duk _ dv^\ _ dum duj dum duk 



dXrn ^ dxi dxk ' dxk dxm dXi dXn 
dum , dum dut dum duk 

\k - -T, 7^ 1- 



ell 



(2.13) 



C^Xfji OXk CfXfji OXi OXfji 

The proof is complete. □ 

lemS Lemma 2.6. Letu{x,t) andp{x,t) be C°° (M^ x [0, oo)) functions satisfying (1.1) 
and (1.2) with fi{x,t) being identically zero for 1 < i < 3, and let {i,k,m) be a 
permutation of (1, 2, 3). The following relations hold for all x and t > 0; 

dXk dXm 

where hi^k is defined by (2.10) and Au, is the Laplacian ofui in the space variables. 
Proof. Prom (2.10) we have 

dhj^k _^ dhj^rn ^ d^Uj _ d^Uk _^ d^Uj _ d'^Um 

dxk dxm dxk'^ dxidxk dxm^ dxidxm 

The right-hand side of (2.14) can be rewritten, by (1.2), as 

d'^Ui d'^Ui dp fduk dujn\ 
dxk"^ dxm^ ddxi\dxk dx^^ 
d'^Ui d'^Ui dp dui ^ 
dxk^ dxm^ ddxi dxi * 
This completes the proof. □ 

As an example of (2.11) let us find another solution to (1.1)-(1.3). 
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lem4 Lemma 2.7. Let bj,aj € M satisfy J2'^=i ^ fii^^t) be chosen identically 

zero for 1 < i < 3. Then 

u° = bism(y^), l<i<3, 

is a smooth, periodic, bounded and divergence-free initial vector field. Let g -.R^M. 
be a smooth function with g{0) = g'{0) = 0. The following family of functions u 
and p satisfy (1.1)-(1.3); 



Ui{x, t) = he^* sin f V ^ + g{t)) - go{t), 

3 

p{x,t) = g'o{t)'^Xj, 



(2.15) |el4 



where g'{t) = go{t) X)^^! ^ and jS = -v X;|=i ^■ 
Proof. Let us write z = X]s=i + foi' brevity. Then 

^ = 0bie^' sm{z) - g'oit) + g'^e^' cos{z), 



3 1 

-vAui = vbiC^* > sin(^), 



= bic'^'siniz) cos{z) j^-- bi9o{t)e^' cos(z) ^ - 



3 1 

-6i5o(i)e''*cos(z) V — 
^ — ' a,- 



since ^ = 0- Inserting the parts to (1.1) shows that 

dui dui dp 

j = l J 

3 3 

by the conditions on g'{t) and /3 in Lemma 2.7. □ 

The simple reasoning leading to the solutions in Lemmas 2.1 and 2.7 is as follows. 
Looking at (2.11) it seems that the leading terms of 

3 

y^^^ dhj^k dUm ^^ ^ dUm dUj dUm Ouk ^ gj^^ 

. . Oxj dx^ryi c)xf^ dxjrj^ dxi Oxjrj^ 



el5 
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should cancel and leave a reminder g{x, t) that can be obtained from the time 
derivative of hi^k- Then there is a first order differential equation 

- vAhi,k + g{x, t) =0 

ot 

which suggests that the solution is exponential and in order to get periodic initial 
values, trigonometric functions were selected. 

In Lemma 2.7 wo first select Uj = (X]s=i ~) where / is a smooth function 
to be determined. This choice automatically gives 

Edhik „ 



3 



because expanding it shows that it has the multiplicative term ^ which is 

zero by divergence- free condition (1.2). The terms 

dum ^ ^ dum duj dum duk 

OXm OXk (-'^m OXj^ OXm 

dUra dUi , dum dUi \ f du,n Ouk du-m duk \ 



V dXm dXk dXk dxm ) ^ dXi dXi dXm ' 

also cancel automatically for the chosen function m. Another way to cancel the 
terms is used in Lemma 2.1. In Lemma 2.1 we originally set Um = hi^k by which 

dhi^k dum , „ 
Um-K ^ — ni,k = 0. 

The remaining terms in the left side of (2.16) are 

dum / _ dUk \ 
dXi V * dXm ' ' 

dum / _|_ du^\ 

V dXrr,. ) ' 



dxk 

The divergence-free condition (1.2) assuming Um = hi^k yields 

Q^-^duj ^ dp / duk \ d / ^ \ 
■f-f dxj ddxi V * dxm ' dxk V dxm ^ 

The form (2.1) is constructed such that it is divergence-free and the term g{x,t) in 
(2.16) can be obtained from (2.16). The way to obtain it is adding a function of t 
to Xi, 1 < i < 3. The basic solution can be further modified by a function g{t) as 
in (2.1) and (2.15). 

Lemmas 2.8 and 2.9, below, generalize Lemma 2.7. Lemma 2.8 shows that 
any periodic smooth function m° = bjg{^j ^) with ^ = can be continued 
to smooth u for zero external force since it can be expressed by its Fourier series. 
Lemma 2.9 generalizes the solution to a case where there are two difii'erent functions. 

lem5 Lemma 2.8. Let fi{x,t) = and 

oo 3 3 

u°{x) =biY^ (cn Sin ( ^ + d„ cos ( ^ ^)) 

71=1 j = l j = l 
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where J2'j^i ^ = 0- '^^^ following functions solve (1.1)-(1.3); 

cx> 3 3 

u,{x, t)=b,J2 e''"'* (cn Sin ( E ^) + ^" E ^) ) ' 

n=l j — 1 j—1 

p{x,t) = 0, 

where /3 = —v '^7^' 

Proof. Direct calculations show that 

dui „ * r> 

which completes the proof. □ 
lem6 Lemma 2.9. Let fi{x,t) = and 

3 3 

u°{x) =6i,isin( V^-) +6j,2sin ( E^-) 
where V,,- -^^^ = for n = 1, 2 anc? 

6,iE^ = — , 6,.2E — • 



ie< /3„ = —v a^- „ /or n = 1, 2. T/ie following two functions solve (1.1)-(1.3); 
Ui{x, t) = bi,ie^^' sin f E — ) + &i,2e'5^* ^^^^ / ^ 

K.,<)=cos(i:5^)»s(i:|^)e<-^ 



l+/32)t 



Proof. Computing uj shows that the term equals — g^' ' . We mention that 
the conditions in the lemma imply ■ — - — = 0, ■ — ^ = y~!,- —3—- Q 

3. Approaches to general initial values 

Let us first notice that the transform by the fimction g{t) in Lemma 2.1 and 
Lemma 2.7 is not a coordinate transform of {x,t) to (x',t') where x'j = Xj + g{t), 
j = 1, 2, 3, and t' = t. The equation (1.1) is not invariant in this transform and if 



then 



-/.(x',t')+5'WE^^^^=0. 
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The transform that is used in Lemma 2.1 and Lemma 2.7, 

u{x, t) u{x' , t') - g'{t), 

3 

p{x,t)^p{x',t') + g"{t)Y,Xj 

i=i 

keeps the initial values w°(x) and J{x, t) fixed if ^(O) = g'{Q) = 0. It is a transform 
that can be done to any solution of (1.1)-(1.3) but it works only for equation (1.1). 

It is certainly not a generally valid coordinate transform. Such should work with 
any equation. It is not valid to think of the transform as 

t) u{x' , t') — g'{t), 
p{x,t) ^p{x',t'), 
Mx,t)^Mx',t')-g"{t) 

because this changes the previously selected force fi{x,t). In fact, what is done 
in Lemma 2.1 is not a change of the coordinate system. The force is kept at the 
selected value at zero. The coordinates are kept at {x, t) as they originally are. The 
pressure is eliminated from (1.1) and (1.2) as in (2.11). The equation (2.11) has 
several solutions for u{x, t) and we find a family of solutions for the initial values of 
Lemma 2.1 and some solutions cannot be smoothly continued to the whole space- 
time. In Theorem 2.4 we notice that it is possible to select a force that picks up 
any of these solutions. 

The equations (1.1)-(1.3) can be solved in a Taylor series form, though summing 
the Taylor series can be difficult. We write 

oo oc oo 

71=0 n=0 71=0 

Equation (1.1) yields 

(n + l)V'7i+l,i + X ^'"'^'"^^ ~ Vl^ll)n,i +Pn- fn,i = 0. 

771=0 ^ 

These equations can be solved recursively by dividing 

'rmj Jn,i 

into two parts 

Vn,i = Vn,i,l + Vn,i,2, 

where Vn^i^i is divergency-free and i'„,i,2 has no turbulence; i.e., it can be obtained 
from some function g{x,t) as 

dg 

dxi 

Thus, what needs to be solved is a system 

(n -I- l)lp„+i^i - vA'lpn,i + Vn,i,l = 0, 

dpn 
Vn,i,2 = -7{ • 

dxi 
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We can see the non-uniqueness of the solution. The division of Vn,i into the two 
parts is not unique: if Ag = 0, then ^ can be inserted to either Vn,i,i or to Vn,i,2- 
We selected a linear symmetric g in Lemma 2.1 in order to have a nice periodic Ui. 

The following approach is another way of finding u{x,t) in (1.1) for a general 
initial vector field vP{x,t). In some cases the method may yield closed form results 
easier than the Taylor series method. If (1.1), (1-2) and (1.3) arc satisfied, p can be 
derived by integration. Let us assume u{x, t) exists. We separate a multiplicative 
part Y{t) such that u{x,t) = Y{t)X{x,t) where the scaling is Y{0) = 1. If there is 
no nontrivial multiplicative factor Y{t) then let us set Y{t) = 1. If (w^)^^ exists 
(e.g. locally for a local solution), and g{x,t) — {u^)~^{X{x,t)) is smooth then 
g{x, t) can be expanded as a power series J2^o ^^dsix) of t and the series converges 
in some small neighborhood of the origin. Since u{x, 0) = vP{x) we have go{x) = x. 
Let us write x'{x,t) = {x'j^{x,t))i<i<3 where 



oo 




s=l 



Thus Ui{x,t) = Y{t)u'i{x'). Let 




du'i dul 
dxk dxi 



and 




Clearly 





and thus 



oo 3 



dgs,r{x) du°{x') 
dxk dx'^ 



dgs,rix) dul{x') \ 
dxi dx'^ J ' 



hiA^)=Y{t)hl,{^')+Y{t)J2t'Y.{ 



s—l r—1 



The interesting term is (here Y' = dY/dt) 





3 
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because it lowers powers of t and allows recursion. Computing the terms in (2.11) 
shows that for some functions Qs,i,k{x,t), it holds 

3 

Q_ dh^^k ^ y^^ dhj^k i9um ^ ^ dum duj dum duk 

dt dXj ^'^ dXm ^'^ dxk dXm dXi dXm 

5=1 r=l ' ' 5=0 

where Qo,i,kix,t) — y(t)/i°j,(a;') + F(t)/Q - ^.(x'). Comparing the coefficients of the 
powers of t individually, and then inserting t = to the equation of each coefficient 
yields, for s > 1, 

3 



dg.A^) dunx) dg,A^) duUx) ^ _ i . ^(^^ q). (3.1) 



y( 

^ V dxk dxr dxi dXj. ) s 

r=l 



From (1.2) we obtain 



=1 



3 



2=1 r=l ' 

= n*)i:^+n*)£<-Ei:^^ 

i=l ' 5=1 i=l r=l ' 

The first term in the right-hand side vanishes because is divergence-free. Again, 
comparing the coefficients of the powers of t individually, and then inserting t = 
to the equation of each coefficient yields, for s > 1, 

= 0. (3.2) 

It seems that solving (3.1) and (3.2) for s = 1 gives gi{x), then we can derive for 
s = 2 and obtain 52(2;) and so on, and that function Qs,i,k contains only terms 
gs',j, s' < s. However, (3.1) and (3.2) do not necessarily determine even gi and this 
approach must be modified. This may be a direction for research how to obtain 
linear systems, like (3.1) and (3.2), from which to continue, but we will not study 
this method more in this short paper. Notice that this approach cannot show that 
a solution u{x, t) exists. 

4. Conclusions 

Theorem 2.2 proves that the solutions to the 3-dimensional Navier-Stokes equa- 
tions for incompressible fluid are not always unique for the initial values and for 
the periodic solutions discussed in Statement D of [1]. This is not surprising, as 
the proof of uniqueness requires periodicity of p{x, t). Periodicity is not required in 
Theorem 2.2 or in [1]. Another (different) counterexample to uniqueness of (1.1) is 
given in [3]. 

Some changes are needed to the official problem setting [1] of the Clay Mathe- 
matics Institute's Navier-Stokes problem. 
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In [1] it is stated that we know for a long time that the initial data u'^{x) can 
be continued uniquely to some finite time. This statement needs a clarification. 
In [2] there is one proof of uniqueness. Temam implicitly assumes that p{x, t) is 
periodic, which [1] does not do. All other easy ways to solve the problem have been 
explicitly excluded. It does not seem to be intended that the pressure p{x, t) should 
be implicitly assumed as periodic. 

If periodicity is not demanded for pressure there are many solutions for zero 
external force, and we can construct a solution that cannot be continued to the 
whole space, as is done in Theorem 2.3. If feedback forces are not excluded, we 
can select the solution from Theorem 2.3 by using a feedback force, as in Theorem 
2.4. Unless Theorem 2.4 is accepted as a proof of Statement D in [1], the official 
problem statement for the millennium problem must be corrected. 

If unique solutions are required then conditions must be imposed on the pressure. 
In the space-periodic case the condition is periodicity of p{x,t), while in the non- 
periodic case growth conditions must be set for p{x,t). It is not enough to make 
these changes only, because if feedback forces are allowed and we require the solution 
to satisfy conditions on p{x, t) then a feedback force can select a solution, which 
does not satisfy these conditions. Such a solution is then a counterexample since it 
is a solution and violates the required conditions. 

The problem statement [1] has an expression a given, externally applied force 
(e.g. gravitation). This is not clear enough to rule out feedback forces since a 
feedback control force is an external force to the controlled system. If the inten- 
tion has been to say that the external force must be similar to gravitation then a 
clarification is needed. Gravitation is time-independent while f{x,t) can be time- 
dependent, thus the similarity to gravitation is not full and the extent of similarity 
needs to be stated. The word given does not clarify the issue since this word has 
many meanings. It could be understood as a weak indication that the external force 
should be selected as a point function and any solution {u,p) is looked for. How- 
ever, the most common meaning of given in this mathematical context is chosen. 
We have to take this meaning to given since otherwise there is a contradiction in 
[1]: the official problem setting claims that the solutions can be uniquely continued 
to some finite time. The only way to understand the uniqueness claim to be correct 
without a requirement of periodicity oip{x, t) is to use a feedback force in the case 
of Lemma 2.1. Thus, some expressions in the official problem setting are unclear 
and understanding given as chosen i.e., not excluding feedback forces, seems to be 
the only one where [1] does not make incorrect statements. The easiest approach 
for improving clarity is to exclude feedback forces. However, while correcting [1] 
one should carefully investigate if there are other unnoticed transforms or other 
problems. It is not a small straightforward correction. 

Let us mention that in a weak formulation of (1.1)-(1.3) it is essential that the 
test function has compact support in M"* x [0, oo) because the presented solutions 
{u,p) with non-periodic pressure cannot be integrated over the whole x [0, oo). 

Acknowledgments. I want to thank the following people: Hcnryka Jormakka and 
Jouni Luukkainen for reading most of the paper and suggesting valuable improve- 
ments to the text, Juha Honkonen and Eero Hyry for useful discussions concerning 
the beginning of the paper, Matti Lehtinen for interesting discussions on many is- 
sues, and the anonymous referees for their comments. Without your help this paper 
would not have been checked. 



14 



J. JORMAKKA 



EJDE-2010/?? 



References 

[1] Charles L. Fefferman; Existence and smoothness of the Navier-Stokes equation, available on- 
line at www.claymath.org. 

[2] Roger Temam; Navier-Stokes equations and nonlinear functions, Regional Conference Series 

in Applied Mathematics 41, SIAM 1983. 
[3] Olga Ladyzhcnskaya; Izvestia Akademii Nauk, SSSR, 33 (1969), pp. 240-247. (in Russian) 

JORMA JORMAKKA 

National Defence University, P.O. Box 07 00861 Helsinki, Finland 
E-mail address: jorma.jormakkafimil.fi 



